Exchange Frequencies in the 2d Wigner crystal 
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Using Path Integral Monte Carlo we have calculated exchange frequencies as electrons undergo 
ring exchanges in a "clean" 2d Wigner crystal as a function of density. The results show agreement 
with WKB calculations at very low density, but show a more rapid increase with density near 
melting. Remarkably, the exchange Hamiltonian closely resembles the measured exchanges in 2d 
3 He. Using the resulting multi-spin exchange model we find the spin Hamiltonian for r s < 175 ±10 is 
a frustrated antiferromagnetic; its likely ground state is a spin liquid. For lower density the ground 
state will be ferromagnetic. 



PACS Numbers: 73.20Dx, 75.10-b, 67.8Jd 

The uniform system of electrons is one of the basic 
models of condensed matter physics. In this paper, we 
report on the first exact calculations of the spin Hamil- 
tonian in the low density 2 dimensional Wigner crystal 
(2dWC) near melting. This system is realized experimen- 
tally with electrons confined at an semiconductor MOS- 
FET's and heterostructures and for electrons on the 
surface of liquid helium Q . 

A homogeneous charged system is characterized by two 
parameters: the density given in terms of r s — a/ao — 
(to* / me){-Ka^p)^ 1 / 2 and the energy in effective Rydbergs 
Ry* = (to* /m e e 2 )Ry where m* is the effective mass and 
e the dielectric constant. Figure 1 summarizes the 2d 
phase diagram. At low density (large r s ) the potential 
energy dominates over the kinetic energy and the sys- 
tem forms a perfect triangular lattice, the Wigner crystal 
H . Tanatar and Ceperley [Q determined that melting at 
zero temperature occurs at r s ~ 37 ± 5. Recent calcula- 
tions H have shown that the low temperature phase is 
free of point defects for densities with r s > 50 but de- 
fects may be present very near melting. At densities for 
r s > 100 the melting is classical, and occurs for temper- 
atures T melt = 2Ry/(T c r s ) where r c « 137 §. 

We determine the spin-spin interaction in the Wigner 
crystal, using Thouless' (§] theory of exchange. Accord- 
ing to this theory, in the absence of point defects, at low 
temperatures the spins will be governed by a Hamiltonian 
of the form: 



= -J2(-i) p jp-Ps 



(i) 



where the sum is over all cyclic (ring) exchanges de- 
scribed by a cyclic permutation P, Jp is its exchange fre- 
quency and Vspin is the corresponding spin exchange op- 
erator. Path Integral Monte Carlo (PIMC) as suggested 
by Thouless || and Roger |Jhas proved to be the only 
reliable way to calculate these parameters. The theory 
and computational method have been tested thoroughly 
on the magnetic properties of bulk helium obtaining ex- 
cellent agreement with measured properties flG]. Rather 



surprisingly, it has been found |Tl| that in both 2d and 3d 
solid 3 He, exchanges of 2, 3 and 4 particles have roughly 
the same order of magnitude and must all be taken into 
account. This is known as the multiple spin exchange 
model(MSE). 
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FIG. 1. Phase diagram. The estimated melting line is 
based on Lindemann's criteria [jjj. The long-dashed line rep- 
resents the Debye temperature. The dotted line is 9, the short 
dashed line is J c . The vertical line is the estimated zero tem- 
perature ferromagnetic (F) transition. 

A WKB calculation of the exchange frequencies in the 
2dWC by Roger || predicted that the three electron J3 
nearest neighbor exchange would dominate, leading to 
a ferromagnetic (F) ground state. It has been recently 
pointed out p2] that if the 2dWC can be stabilized by 
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disorder to higher densities, a transition to an antifer- 
romagnetic(AF) phase will occur. Although the 1/r in- 
teraction is characterized as 'soft', in the 2dWC the two 
particle pair correlation function for the electron system 
is quite similar to that of solid 3 He supporting the idea 
that multiple exchanges could be important in the 2dWC 
near melting. 

The method of determining the exchange frequency in 
quantum crystals has been given earlier One 
computes the free energy necessary to make an exchange 
beginning with one arrangement of particles to lattice 
sites Z and ending on a permuted arrangement PZ: 

F P (tf=Q(P,0)/Q(I,0)=taah(Jp(p-po)). (2) 

Here Q(P,[3) is the partition function corresponding to 
an exchange P at a temperature I is the identity 

permutation. Note that these paths are of "distinguish- 
able" particles since Fermi statistics are implemented 
through the spin Hamilitonian in Eq. (1). We deter- 
mine the function Fp(f3) using a method which directly 
calculates free energy differences and thereby determines 
Jp and /3o- The only new feature with respect to the 
3 He calculations is the method of treating the long-range 
potential. We use the standard Ewald breakup |J and 
treat the short-range part using the exact pair action 
fiof and the long-range k-space term using the primi- 
tive approximation. We used a hexagonal unit cell with 
periodic boundary conditions, most calculations contain- 
ing 36 electrons. Checks with up to 144 electrons did 
not change the results within the statistical errors. The 
number of particles is not as important in the 2dWC as 
in solid 3 He because the 1/r interaction suppresses the 
long wavelength charge fluctuations. 



tion functions would be signaled by Fp{j3) not described 
by Eq. (2). Except for r s < 50 where our calculations are 
too unstable to make definite predictions, we observed no 
problems of convergence. 

Figure 2 shows the ring exchanges considered here. Ex- 
cept near melting these exchanges give rise to most of the 
thermodynamic properties. Note that we consider the 6- 
particle parallelogram (6p) exchange, which is not taken 
into account in solid 3 He. We have also calculated sev- 
eral 2-5 particle exchanges having next-nearest neighbor 
exchanges and all possible 6 particle nearest-neighbor ex- 
changes for 50 < r s < 75, but because their magnitudes 
are much smaller, we do not report those results. 




2 3 4 5 6h 6p 

FIG. 2. Considered exchanges. 

We performed calculations for the densities r s = 
{45, 50, 60, 75, 100, 140, 200}. We found accurate results 
using a "time-step" for the descretized imaginary time 

3/2 

path integrals of r < 0.3r s and extrapolated to the 
t = limit using J p (r) = J P (0) + JpT 3 . The inverse 
temperature (3 in Eq. (2 )must be larger than the values 
determined by melting and twice the exchange "time" 

3/2 

(3q ~ 5ry . We have determined the exchange frequen- 
cies with an accuracy between 1.5% and 6% independent 
of their magnitude or the number of exchanging electrons 
though the computer time increases with the number of 
exchanging electrons. Breakdown of Thouless' theory 
caused by many states contributing to the ratio of parti- 



1/2 

FIG. 3. Exchange frequencies versus r s . The solid line 
is 10 -3 of the kinetic energy. 

Our calculated exchange frequencies vary rapidly with 
density as shown in fig. (3). One can see that they are 
much less than the zero point energy of the electrons, thus 
justifying the use of Thouless' theory. The WKB method 
]ll"| where one approximates the Path Integral in Eq. (2) 
by the single most probable path, explains most of this 
density dependence. In the 2dWC, the WKB expression 
for the exchange frequency J9|,[l2| is: 

J P =A P (r s )b 1 J 2 r-^e- b ^ /2 . (3) 

1/2 

Here bpr s is the minimum value of the action inte- 
gral along the path connecting PZ with Z. The 3 parti- 
cle exchange exponent is the smallest indicating that as 
r s — > oo, J3 will dominate and the system will have a 



2 



ferromagnetic ground state. However, note that in fig. 
(3) J2 > J3 for r s < 90. We observe a more rapid in- 
crease in all the exchange frequencies than predicted by 
Eq. (3) for r s < 100 caused by non-linear fluctuations 
about the classical path. Near melting the spread in ex- 
change frequencies is much less than in the low density 
limit. 
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FIG. 4. The prefactor of the exchange frequencies, Ap , 
see Eq. (3) versus density. All prefactors go to a constant 
for r 3 > 100 indicating that exchange frequencies are domi- 
nated by a single path (WKB limit). The curves are fits as 
described in the text. 

Figure (4) shows density dependence of the prefactor, 
Ap. The prefactors are larger than unity as expected 
|l2| (our definition of A p differs from ref. jL2| by a factor 
of 1.63/-\/27r). We observe that the prefactor is roughly 
constant for r s > 100 but rises rapidly as the system 
approaches the melting density. We are unaware of any 
work that definitively establishes the form of the correc- 
tion to the WKB formula. A good fit was obtained with 
the function Ap(r s ) — Ap(0)(l + (rp/r s ) 3 ) to determine 
bp, Ap(0) and rp. The exponents are close to the WKB 
values, only differing significantly for the 2-particle ex- 
change. 

The prefactor for J2 is twice as large as the other 
exchanges. One might think that since one cannot ob- 
serve the difference in pair exchanges in the initial and 
final spins, the sense of the exchange would not mat- 
ter. Our results show that the 2 orientations of 2 par- 
ticle exchange, clockwise and counter-clockwise, should 



be counted separately. A recent calculation Q directly 
evaluates J2 in the high density regime, 5 < r s < 45, by 
numerically solving for the difference between the even 
and odd parity 2-electron energies. In those calculations, 
the spectator electrons were fixed at their lattice sites 
and the 2 exchanging electrons had vanishing wavefunc- 
tions outside a rectangle centered around the exchange. 
At r s = 45 a direct comparison shows their result is 2.3 
times larger than that obtained with PIMC. This result 
is expected because of the additional localization caused 
by the spectator electrons fluctuating into the exchanging 
region. 

Having determined the exchange frequency, one is left 
with the spin Hamiltonian of Eq. (1). This is a non- 
trivial many-body problem which we will not discuss 
in detail here. For spin 1/2 systems, J2 and J3 con- 
tribute only with a nearest neighbor Heisenberg term: 
J| ff = J2— 2J3. This term is negative (ferromagnetic) but 
approaches zero near melting. For convenience we use J4 
as a reference to fix the overall scale of the magnetic en- 
ergy. Neglecting Jq p , the Hamiltonian has 3 remaining 
parameters Jf /J4, J5/J4 and Jqu/Ja- The dependence 
of these ratios on density is shown in fig. (5). 




-5 -4 -3 -2 -1 

JfA 



FIG. 5. Spin Phase diagram as a function of exchange ra- 
tios. The dotted line is the flow of spin hamiltonian space 
versus r s (lower numbers); also shown the are estimated 
values of Jen/Ji (upper numbers). The solid lines are the 
limit of the ferromagnetic phase according to ED |l7j at 
Jeh/Ji = {0.2,0.4,0.6}. The 2dWC crosses into the F re- 
gion for r s ~ 175. The (*) are empirical estimates of the spin 
Hamiltonian of 2d 3 He at several densities fllSJ. 
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High temperature series expansions |lq| determine the 
the specific heat GV and magnetic susceptibility Xo/x f° r 
temperatures ksT ^> Jp. The susceptibility is given by: 
Xo lx = T - 9 + B/T . . . and the specific heat C v /Nk B = 
(3J C /2T) 2 + . . . where the Curie- Weiss constant is given 
by 9 = -3(J 2 0ff + 3J 4 - 5J 5 + 5/8J 6h + 15/4J 6p ) with a 
quadratic expression of the J's for J c . These two con- 
stants which set the scale of the temperature where ex- 
change is important are shown as a dotted and dashed 
lines on fig. (1). Both 9 and J c decrease very rapidly at 
low density showing that experiments must be done at 
r s < 100 if spin effects are to be at a reasonable temper- 
ature, e.g. T c > O.lmK (assuming e = m* = 1). Note 
that 9 changes from positive to negative at r s w 130. 

The zero temperature state can be studied by exact 
diagonalization (ED) of an N site system. (The present 
limitation is N < 36.) Here we focus on the zero tem- 
perature transition from a ferromagnetic to an antiferro- 
magnetic state. An estimate |l7| of the critical transition 
parameters is shown on fig. (5). The ferromagnetic phase 
is obtained only at very low density: we estimate the the 
F-AF transition for the 2dWC will occur at r s = 175±10. 
(Note that this estimate does not include the effect of Jq p ; 
this will increase the stability of the antiferromagnetic re- 
gion to roughly r s = 200.) 

At higher density, the frustration between large cycle 
exchanges (4-6 body loops) results in a disordered spin 
state. For example, the point (J| ff /J4 = —2, J 5 = 0, 
J6h = 0), close to the parameters at r s = 100, is a spin 
liquid with a gap to all excitations |T^] : The spin liquid 
properties can be understood from a resonance valance 
bond model with ordering into spin-1 diamond plaque- 
ttes. The breakup of the plaquettes is responsible for 
a low temperature peak on the specific heat. A second 
high temperature peak develops at a temperature T ~ J c , 
shown as a dotted line on fig. (1). 

We find that at higher densities (r s < 100), the trajec- 
tory of the MSE models parallels the F-AF phase line, 
with the possibility of a re-entrant ferromagnetic phase 
for r s < 40. Note that QMC calculations (T|] of the nor- 
mal fermi liquid at r s = 30 show that the ferromagnetic 
phase has a slightly lower energy that the unpolarized 
phase. Hence both the high density 2dWC and the low 
density electron fluid are characterized by a spin Hamil- 
tonian which is nearly ferromagnetic. 

We note a remarkable similarity between the exchange 
parameters of the 2dWC to those extrapolated from mea- 



surements of the second layer of 3 He absorbed on grafoil 
|]l9f as shown in fig. (5) . The existence of these two re- 
lated spin systems should allow a fuller understanding of 
this remarkable spin liquid. Such an equivalence could 
arise from an underlying virtual vacancy-interstitial (VI) 
mechanism [ fLl| giving rise to ring exchanges. In this 
model the prefactor of the exchange is controlled by the 
rate of VI formation (non-universal) but the ratios of the 
various ring exchanges arises from geometry of the tri- 
angular lattice and from the attraction of point defects 
(universal). We have recently determined || that the VI 
formation energy vanishes at melting in the 2dWC. This 
is consistent with the fact that the various Jp increase 
rapidly near melting. 

In summary, we find that the magnetic interactions are 
characterized by a frustrated spin order. Application of a 
magnetic field pcfl transforms the exchange frequencies to 
j pe 27TteB ± a p/ h where a p is the area of the exchange (see 
table I) and B± the magnetic field. Experiments with 
magnetic fields will allow exploration of this Aharonov- 
Bohm effect and thereby provide direct information on 
ring exchanges. 

The semiconductor realizations of the 2dWC have sig- 
nificant disorder which can stabilize localized electronic 
states at higher densities than in the clean system pj[ . 
One can also stabilize the Wigner crystal at higher den- 
sities using bilayers p2| . Exchange frequencies in those 
systems including effects of layer thickness and the ex- 
change properties of point defects could be calculated 
with the present method. 

This research was funded by NSF DMR 98-02373, the 
CNRS-University of Illinois exchange agreement, support 
by Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo (FAPESP) and the Dept. of Physics at the Uni- 
versity of Illinois. We used the computational facilities 
at the NCSA. 
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1.612 


2.11 


44 
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1.525 


0.91 


47 


2.07(1) 


4 


1.656 


1.42 


45 


3.07(2) 


5 


1.912 


1.39 


60 


4.11(2) 


6h 


1.790 


1.12 


59 


7.04(3) 


6p 


2.136 


6.24 


38 


5.09(2) 



TABLE I. Fits for the exchange frequencies. The path 
area, ap , was calculated at r 3 = 60 and it is in units of trian- 
gle areas. 
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